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ABSTRACT 



The ability to measure and predict atmospheric 
turbulence affecting laser beam propagation is a major 
concern when considering military applications. Such a 
method using a telescope, high resolution television camera, 
video tape recorder, digital storage oscilloscope, and 
calculator system has been devised, tested and utilized. A 
laser beam signal is recorded on video tape for further 
processing. This signal is displayed, stored and digitized 
using a Tektronix 453 digital storage oscilloscope. The 
digitized signal is sent to a Hewlett-Packard 9825 comparing 
system for Fourier transform aaalysis and determination of 
the refractive index structure constant, C n 2 • Several 
trials were conducted using He-Ne and G a - A s lasers. The 
results demonstrated good correlation with theoretical 
predictions as well as previously analyzed data. 
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I. 



INTRODUCTION 



A. BACKGROUND 

The increasing use of lasers and laser technology for 
military applications has brougat about a need for analysis 
of the laser beam in its environment, the turbulent 
atmosphere. A pro jet t at Naval Postgraduate School dealing 
with this subject is continuing and is the main topic of 
this thesis. 

The patterns produced by lasers on targets have inherent 
problems that include broadening, beam wander and intensity 
fluctuations brought about by turbulence in the atmosphere. 
These effects of atmospheric turbulence on laser propagation 
have been wall determined [Ref. 1]. In terms of the Fried 
model, [Ref. 2], C n 2 , the refractive index structure 
constant, has been determined to adequately express the 
Modulation Transfer Function (MTF) , or the Mutual Coherence 
Function (MCF) for the atmosphere. 

A system has been developed that provides a measurement 
of C n 2 for atmospheric turbulence along the optical path 
through which a laser is propagated. The system employs a 
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vidicon and telescops as the detector and a distant laser as 
the source. This system duplicates the slit scanning system 
presently in use at Naval Postgraduate School [Ref. 2]. 
Tests using the vidicon equipment have been previously 
completed with values for C n 2 on the order of 10 -15 m -2 / 3 
being obtained [Ref. 3]. Measurements have been made using 
the same experimental set-up, except that the Tektronix 468 
Digital Storage Oscilloscope are used for digitization 
instead of the Quantex DS-30 Digital Video Processor. These 
measurements have demonstrated values of 3 n 2 of comparable 
accuracy. 

B. PROBLEM 

3y using the digital storage capability of the Tektronix 
468, the previously recorded laser signal is used as an 
input and is evaluated by modifying the program developed by 
Craaer [Ref. 4]. h brief overview of the procedure is 
described below. Detailed explantions of the experiment and 
computer program are contained in Chapters III and IV, 
respectively. 

The approach taken is basically the same as that of 
Crager, but because of the unavailability of a disk ROM it 
is necessary to store both raw and processed data on 
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magnetic taps. Ths data are loaded and recorded as 

necessary by the HP 3325. The basic assumption permitting 
analysis is that the horizontal TV scan line through the 
laser spot is considered to accurately mirror a point spread 
function of the image. Final analysis using the 463 



oscilloscope has shown 


this to 


be 


a valid hypothesis. 


The 


digitized 


data from 


the 468 


D sci 


iioscope agree with 


the 


previously 


measured ia 


ta which 


i S3 


many pixels of televi 


sion 



data digitized by the DS-30 to express the point spread 
function . 

The sequence of analysis is that a video tape recording 
is made of the TV image of the propagated laser beam. The 
output of the video recorder is sent to the Tektronix 468, 
where the derived TV scan line is digitized, averaged, and 
stored. The HP 9825 records the digitized data and produces 
a line spread function (LSF) , by integrating the point 
spread function. 



L3 F (X) 



y (max) 
J?SF(r) dy 
y (min) 



( 1 . 1 ) 



where 



r = V< x 2 + y 2 ) ' 

and computes the Fourier transform of the LSF. The 
diffraction limited Fourier cransform of the optics is now 
computed if the MTF of the optics has not been previously 
measured experimentally. Next, the program finds the MTF of 
the atmosphere by dividing the Fourier transform of the LSF 
by the Fourier crar.sforn of the optics. Finally, by curve 
fitting, the program computes a single value for C n 2 . 

The program now predicts the size of a laser spot on a 
target using the calculated value of C n 2 . This value of C n 2 
is used to calculate an HTF of the atmosphere which is chen 
multiplied by the Fourier transform of a source, and one 
Fourier transform of the optics. The program then 
calculates the inverse Fourier transform of the produces of 
the above and uses the Abel transform to give the angular 
point spread intensity distrioution. From this data, ohe 
fraction of energy as a function of the total energy within 
a given radius R is calculated. 
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II. 



THEORETICAL -DMSIDERATIDNS 



The theory of laser beam propagation through a turbulent 
medium has been explained by Crittenden, and others, and is 
re-emphasized here for continuity purposes [Ref. 2]. Since 
the effects of turbulence on laser weapons is of major 
concern, measurement and prediction capability for these 
effects on laser beams is entirely relevant. 

Due to the existence of reciprocity, these ideas apply 
to either laser designators or imagers [Ref. 5]. In the 
study of the theoretical modeL by Fried, the effects of 



atmospheric 


turbul en c e 


are investigated [Ref. 6]. 


This 


model uses 


the idea of a 


long term optical transfer 


f unc tier. 


(OTF) when 


c onsideri ng 


at mo son eric turbulence. 


The long 



term OTF results from simply taring an image of sufficiently 
long term which sees effectively all possible turbulence 
conf igurati cns. 

When observing tie effect of diffraction as the result 
of the finite aperture of a point source, it can be seen 
that the resulting image is not uniform. Considering the 
point spread function to have the same shape in the image 
plane regardless of its position, the image function is rhe 
result of the convolution of the source function and one 
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optics diffraction function. The convolution theorem, as 
described in Fourier transforn theory, yields 

i (v (x) , v (y )) =H (v(x) , v (y) | *o (v (x| , v (y) ) *M (v (x) , v (y) ) (2.1) 

where 



i(v(x),v(y)) = Fouriec transforn of the image function 



H {v (x) ,v (y) ) = 



cun er 
: unct ion 



C3 



( ^ r*N ^ A - G • 



o(v(x),v(y)) = Fourier transforn of the object function 

M (v (x) , (v (y ) ) = Modulation Transfer Function of 

the atmosphere 

v = spatial frequency 

A point source such as a Laser car. be analyzed in two 
dimensions using a point spreal function. Through Fourier 
transform theory, the image point spread function is 
transformed into a two dimensional OTF of the optical 
sianal. This problem may be simplified by scanning the 
image point spread function using a vidicon or slit-scanning 
system. In applying the convolution theorem, the Fourier 
transform of the image point spread function is multiplied 
with the Fourier transform of the optical system resulting 
in the Fourier transform of the overall system. The Abel 
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■transfer described by Griem, is applies to this r-sti : to 
re-transf or a the on e -dim ensio a al image LSF into a two 
dimensional image PSF [Ref. 7], 

As demonstrated by Crittenden, and others, a numerical 
value for C n 2 may be obtained by curve fitting using the 
following: 

M = exp(-21.49*C n 2 *Z*f 5 / 3 ^A~ 1 / 3 ) (2.2) 

where 

M = MTF of the atmosphere 
Z = Range in meters 

f = F* v = Angular spatial frequency in cycles/radian 
F = focal length of the optical system 
v = linear spatial frequency in cycles/meter 
A = wavelength in meters 

C n 2 is obtained by a linear regression of 

In (M) = -21.49*C n 2 *Z*f s/3*^-x / 3 (2.3) 

where C n 2 is the only parameter. 
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III. 



EXPERIMENTAL PROCEDURE 



The experiment is performs! using two different laser 
sources, helium-neon (He-Ne) and gallium- ar sen ids (Ga-As). 
Figure 3.1 shows a block diagram of the experimental set-up 
used for both lasers. 

The measurement of C n 2 along the optical path is made by 
using a viaicon and telescope at the far end of she corridor 
in the basement of Spanagel Hall. Atmospheric turbulence is 
produced by nine overhead hot air ventilators. The optical 
equipment similar to that described in reference 3, consists 
of a 6 inch diameter Cassegrain telescope with a 90 inch 
focal length [Ref. 3]. 
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Fig. 3. 1 Laser Experimental Set-Up 






analysis is accomplished using a Tektronix 468 Digital 
Storage Oscilloscope and a Hewlett Packard 9825 Comparing 
System. The linearity of the Panasonic tape recorder was 
demonstrated in Crager's thesis. This is also confirmed by 
observing the images as real time and recorded 
displays [3ef. 3], 

The 468 oscilloscope is set-up in the following manner. 
An input signal to channel A is located on the scope using 
the non-storage mode. The signal is then displayed using 

the A INTENSITY switch of the horizontal display. This 
intensified zone is used to position the 3 sweep (delayed) 
to the desired location within the A sweep interval to 
obtain an expanded view of a waveform for examination. Once 
the waveform is centered on the scope, the horizontal 
display is switched to 3 DLY'D. This is done to facilitate 
the digital storage circuitry time base by using the setting 
of the B TIME/D IV switch. 

The waveform is expanded in time by decreasing the A 
TINE/DIV switch setting and moving the waveform back to the 
middle of the scope using the delay time position control 
dial. When a representative waveform is obtained, the B 
TIME/DIV switch is then used in conjunction with the delay 
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time position control dial to expand the tins scale until 
the single central maximum of the waveform is centered ana 
one horizontal sweep is displayed. 

Once this wavefocn is satisfactorily obtained, the AVG 
storage mode is selected. The 468 will average the input 



signal for a selected number of 


sweeps and display the 


accumulated waveform. All data 


for this thesis use 32 


sweeps for each average. The 468 


is now ready to transfer 



data when interrogated by the HP 9825. 

The computer is the controller for all interfacing 
operations with only a minimal amount of operator 
interaction. The operator interface is mainly to ensure 



that the equipment is properly set 


up and to select if data 


are to be plotted. Digitization 


of the analog signal is 


accomplished by the Tektronix 45"?. 


Processed waveform data 



are transferred from the microprocessor memory to the 
Storage Display RAM [Ref. 8]. 

The controlling program of the HP 9825 interrogates the 
468 via the IEEE 488 interface bus. The 468 receives the 
data request from the controller and sends the waveform 
message, both preamble and data. The waveform message is 
stored in the calculator memory for further processing. 
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When the message is completed, the 468 concludes with an er.d 
of instruction terminator and the controller takes control 
of the bus again. 

Data processing of the digitized signal begins with the 
data being stored on magnetic tape for further use by the 
computer. The data are stored, processed and plotted by the 
main program. Subroutines are called as necessary for their 
specific uses. The two signal waveforms are recorded and 
processed. First, the signal from a laser beam incident at 
the aperture of the telescope is recorded. Figures 3.2 and 
3.3 show these recorded data waich are referred to as the 
source data. Second, for calibration, the signal from a 
laser beam with a diffraction grating in place at the 
aperture of the telescope is recorded. The entire system is 
calibrated by using a grating that consists of closely 
spaced vertical bars in front of the telescope. This 
produces a diffraction pattern in the image plane and allows 
for the calculation of the scale factor. These recorded 
data are referred to as the scale data and are plotted in 
Figures 3.4 and 3.5. 

During the recording process, the signal from the laser 
beam is recorded on video tape while the signal is observed 
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Fig. 3.2 Plot of Source Data for He-Ne Laser 



N 




AX 1 SN3XN I <33Z 1 1UWH0N 



22 



Fig. 3. 3 Plot of Source Data for Ga~As Laser 
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Fig. 3.4 Plot of Scale Data for He-Ne Laaer 
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Fig. 3.5 Plot of Scale Data for Ga-As Laser 



on nhe U58 osciiicsoo pe. This procedure is important in 
data taking since r a s e r axrgnmsni anu texescops roc u sang are 
often very sensitive to minor movements of either. Upon 
completion of recording the data on video -ape, the 
experiment is concerned with the compilation and analysis of 
data. The first item to be determined is the scale factor. 
Once the scale factor has been determined, it will remain 
constant throughout the calculi- ions, unless tne telescope 
focal length is changed. 

The scale factor is calculated in the following manner. 



The number 


of points between the 


peaks in the plot 


of the 


scale data 


(diffraction grating 


in place in front 


of the 


telescope a 


perture) is measured. 


When the distance 


bet ween 


peaks, the 


spacing between the bar 


s in the grating. 


and the 


wavelength 


cf the laser are known 


, the scale factor 


ca n be 


calculated 


from the relation 







sin (theta) = X /d 


(3. 1) 


sin (theta) ~ theta 


(3.2) 


(small angle approximation) 




sf = (theta) /ndp 


(3.3) 


sf = A / (d*ndp ) 


(3.4) 
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wn er e 



sf = scale factor in radians per point 
A. = wavelength in meters 

a = spacing between lines on grating in tneters 

ndp = number of points between central maximum 
and first order diffraction peak 

This standardizes the data for the abscissa of the plots in 
radians per point. 

The program now takes the transferred data and computes 
a point spread function as shown in Figures 3.5 and 3.7. 
Figures 3.3 and 3.9 represent the point spread function 
after integration using equation (1.1) to obtain a line 
spread function. Next, the Fourier transform of the line 
spread function is calculated. These curves are plotted in 
Figures 3.10 and 3.11. 

The diffraction of the optics is computed and plotted as 
shown in Figures 3.12 and 3.13. Following the same method 
as before, the line spread function of the optics function 
is calculated and -plotted in Figures 3.1'u and 3.15. The 
program now takes the Fourier transform of the optics 
function LSF. These data are plotted in Figures 3.16 and 
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Fig. 3.6 Plot of Point Spread Function for He-Ne Laser 
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Fig. 3.7 Plot of Point Spread Function for Ga-As Laoer 
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The Fourier transforms of the system and optics are 
divided point by point to yield the MTF of the atmosphere, 
as shown. in Figures 3. 18 ani 3.19. The MTF of the 
atmosphere will be used to calculate the long-term value of 

Cn 2 - 

C n 2 is calculated by linear regression of In (MTF) versus 
f s / 3 . The slope of the regression is proportional to C n 2 . 
Two simultaneous equations are solved using Cramer's rule 
and the matrix ROM of the HP 9825. The equations used are 

a*X + b*I = Y (3.5) 

a *X 2 + b*X = Y*X (3.6) 

where 

a = slope of the curve 
b = intercept 

X = summation of (I*sf ) ** ( 5/3) 

X** 2 = summation of (I *sf) ** (1 0/3) 

I = point number (total number = 256) 

Y = summation of natural logarithm of M7F 
Y*X = summation of products of two values 
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C n 2 is obtained from the aoove information using equation 



(2.3). This will yield 

C n 2 = a/(-21. 49*Z*f ** (5/3) * **(-1/3)) (3.7) 

The program now goes to the prediction phase after 
calculating C n 2 for the atmosphere. If a Gaussian 
distribution for an input source and a value of the standard 
deviation are assumed or known, the resulting source 

function can be calculated. A plot of the computed source 
is shown in Figure 3.20. It can be seen from this plot that 
the half- width at half-maximum is approximately 4 

microradians. As before, the line spread function of the 
computed source is calculated and plotted. This is shown in 
Figure 3.21. The Fourier transform of this data is 
calculated and plotted in Figure 3.22. 

The program now multiplies the Fourier transform of the 
computed source with the transform of the system (including 
the atmosphere) and plots the result in Figures 3.23 and 
3.24. Plots of the inverse Fourier transform are shown in 
Figures 3.25 and 3.25. Next, the Abel transform is computed 
and plotted in Figures 3.27 and 3.23. The Abel transform, 
as described in Chapter II, transforms a one-dimensional 
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sp read 



line spread function into a two-dimensional point 
function. Finally, one fraction of power inside a circle of 
radius R is calculated and plotted as shown in Figures 3.29 
and 3.30. This is the fraction of power that one would 
expect to be incident on a target using the measured value 
of atmospheric turbuLence as an input parameter. A detailed 
analysis of the computer program is the suoject of the next 
chapter. 
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Fig. 3. 8 Plot of Line Spread Function for He-Ne Laser 
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Fig. 3.9 Plot of Line Spread Function for Ga-As Laser 
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Fig. 3. 10 Plot of Fourier Transform of LSF for He-Ne Laser 
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Fig. 3.11 Plot of Fourier Transform of LSF for Ga-As Laser 
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Fig. 3. 12 Plot of Optics Diffraction for He-Ne Laser 
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Fig. 3. 13 Plot of Optics Diff root ion for Go -As Laser 
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Fig. 3. 14 Plot of LSF of Optics Function for He-Ne Laser 
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Fig. 3. 15 Plot of LSF of Optics Function for Ga-As Laser 
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Fig. 3. 18 Plot of MTF of Atmosphere for He~Ne Laser 
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Fig. 3. 19 Plot of MTF of Atmosphere for Ga-As Laser 
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Fig. 3.20 Plot of Computed Source 
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Fig. 3.21 Plot of Line Spread Function of Computed Source 
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Fig. 3.22 Plot of Fourier Transform of LSF of Computed Source 
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Fig. 3.25 Plot of Inverse Fourier Transform for He-Ne Laser 
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Fig. 3.26 Plot of Inverse Fourier Transform for Ga-As Laser 
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Fig. 3.27 Plot of Abel Transform for He-No Laser 
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Fig. 3.28 Plot of Abel Transform for Ga~As Laser 
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Fig. 3.29 Plot of Power Fraction Inside Circle of Radius R for He-Ne Laser 
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Fig. 3.30 Plot of Power Fraction Inside Circle of Radius R for Ga-As Laser 



IV. 



COMPUTES PROGRAM IS A 



GORITHMIC LANGUAGE 



The program presented here is written in the algorithmic 
language as outline! by Grahan [Ref. 9]. It is intended 
that the language in use here be machine independent. 
Comments appear in parentheses. 



Alaorithm_ Laser 
(Input which laser is used.) 

Input Laser_Navelength 
ESTRING <== Type_of _Laser 

(Input if the particular laser source data are recorded. 
The following lines determine if the laser source data are 
already recorded. If the source data are not recorded then 
the program goes to the transfer subroutine to allow 
interfacing with the 468 to effect a data transfer. The 
parameters passed to the Plot subroutine have the following 
significance: the first parameter determines the horizontal 

index, for example 512 points or the number of microradians. 
The second parameter tells the computer which file the data 
are recorded on so the data can be loaded into memory. The 
final parameter determines the labeling of the horizontal 
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_ ii 



or 



axis either in data points, microra dians, 
microradian .) 

If Laser_Source_no t_recori9d then 
Call Transfer (Source_Fi le) 

Call Plot (51 2 ,Source_File, 0) 

End If 






(Input if the particular laser scale data are recorded. 
Scale data are recorded with the diffraction grating in 
place. If the scale data need to be recorded, the 
subroutine Transfer is called and the scale data are stored 
on a separate file.) 

If Lase r_S.cale_not _rer orded then 
Call Transfer (Scale_File) 

Call Plot (512 ,Scale_file,0) 

End If 

Input Scaie_Factor 

(Calculation of the scale factor is detailed in Chapter III. 
Comparison of the time base of the recorded data is 
necessary to ensure that the data recorded have the same 
time scale on the oscilloscope. If the time scales are 
different the program stops and displays an error message. 
A file is loaded, the waveform preamble is searched for the 
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“X'JKITS" or "XINCR" oscilloscope settings end the 
particular data are stored. This is done for each of the 3 
recorded files: source, scale, and collimated data. 

Collimated data refer to the same laser recorded under 
controlled conditions. The laser beam is projected in a 
collimator and recorded on tape. This data could be used to 
remove any non-linearities of the video recorder, electronic 
equipment, or receiving optics. Data elements are compared 
to ensure consistency of inputs.) 

Input Laser_source 

3<== po s(Laser_So u rce, • H : 1 ) 

R 1 <== val ( • B+2» ,3+4») 

3<== po s (laser_Sou res, 1 IT: } 

ASTRIN G <== Laser, Source ( • B + 3 1 , 1 3+4 » ) 

Input (L aser_Sca lei 

3<== po s (L a ser_ Sc a le , 1 R : ’ ) 

R2<== val ( • B+2 1 ,3+ 4» ) 

B<== po s(Laser_Soa le, ‘ IT: ' ) 

BSTRIN G<== Laser_3 cale ( 'B+3 • , ' 3+4 •) 

Input (Colli mate d_3 ource) 

B<== po s(Collimate d_So urea , 1 R : ' ) 

R3<== val ( ' B+2' ,» B+4) 

B<== pos (Collin ate d_So urea, 1 IT: ') 
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CSTR iNG <= = Colii:nnfed_Sourc5('B + 3', * 3 4 * ) 

If R 1 = S 2=R 3 then 

If ASTRIN3 = BSTRING = CSTRING then 
Input SO ATA ( File_N umber) 

End If 

Else 

Output 'ERROR IN TIME SCALE' 

Step 
End If 

A<==B<==C<==0 

(Array is initialized to 0) 

Do for I = 1 to 51 2 

YD ATA (I) < = = R DATA (I) 

End Do 

(This section takes the data from the waveform message as 
transferred from the 468 and removes any DC background from 
an image as well as subtracting the minimum value in order 
to "zero" the data. This makes y=Q as the starting point 
for both calculation and plotting purposes. Data up to 51 
and greater than 463 are summed and stored in B. 

This represents the "wings" of the curve. This 
process removes the DC background from an image.) 

Do for I = 1 to 51 2 
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i< = A + YD A IA ( I) 



If I < 51 than 

B <== YDATA(I) +3 
End If 

If I > 463 then 

3<== YDATA(I) + B 
End If 
End ‘Do 

(Output for debugging) 

Ourpuo A, 3 

(Average of the sum of data values in the "wings".) 

D<== B/100 

(This is the total sum of all data less the "wings" value.) 
A<== A- 51 2 *D 

(All data are summed and stored in C. Each data 
is compared with the previous total divided by 2. 

This will find the index value, I, for which 
the maximum value exists.) 

Do for I = 1 to 512 

Y D AT A ( I) <== Y DATA ( I) -D 
C<= YDATA (I) + E 

If C < A/2 then 
J<== I 
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E r. d If 
End Do 

(Output for debugqiig.) 

Output C,J 

(The curve is shifts! so that ths maximum value starts 
at the origin of ths horizontal axis. This is called 
the point spread function and is recorded as RDATA.) 
Do for K = 1 to 256 
I DATA (K) < = = 

YD AT A ( ( J -K + 1) mod512 + 1| +YDATA ((Ki- J-1) mod512*1) 

End Do 

Do for I = 257 to 512 
RDATA (!)< = = D 
End Do 

Do for I = 1 to 256 

RD ATA (I) <== IDATA (I) 

End Do 

(Determines if tabular or plotted data ars desired.) 

If PSF_OUIPUT_waated then 
Call Table (RD ATA| 

End If 

If PSF_ PL3T_wan ted then 

Call Plot (32 , RDATA, 1| 
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End If 



(Calculates the line spread function from the point 
spread function.) 

Call LSF 

(Determines if tabular or plotted data are desired.) 
If LSF_ OUT? UT_ wanted then 
Call Table (RDATA) 

Else If LSF_PLOT_w anted than 
Call Plot (32, RDATA, 1) 

End If 

(Calculation of Fourier Transform of LSF.) 

Input RDATA 

(The maximum value of the array IDATA is at 1=1. Thi 
is put back into array RDATA and folded over so that 
IDATA (257) to IDATA(5 12) is the mirror image of 
IDATA (1) to IDATA (255 ). Then this array is recorded 
as RDATA and the Fourier transform calculated.) 

Do for I = 257 to 512 
RDATA (I) < = = 3 
End Do 

Do for I = 1 to 256 

RD ATA (514-1) < == R DATA (I) 

End Do 
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P.DATA (2 57) < = = 3DAPA{256) 

Call F X FOE M 

(Determines if tabular or plotted data are desired.) 

If FXFORM_OF_LSF_DUTPaT_wantad then 
Call Table (RD AT A) 

Else If FXFORM_OF_LSF_PLOT_wanted then 
Call Plot (32 r RDAT A, 2) 

End If 

(Dia of Cbscur/Dia of Dbj Lens in meters , 3 ATIO= . 064/. 164) 
(This is the measurement of the ratio of 
diameters of obscuration to telescope.) 

(R1= Scale of data in microradians=1. 60) 

(Diameter of telescope Obj Lens in meter s , 0BJECT= 

0.164 me ters) 

If Laser_Source = He-Ne than 
Wa velength<== 6.328*e-7 

Else 

4a velength<== 9.05*e-7 
End If 

(This begins the calculation of the Fourier transform 
of the optics using the Air/ function. All of the 
data regarding the optics are recorded on a separate file.) 
(These next four statements fix the constants in the 



52 



argument of the optics function.) 



R 1 <== 3cale_Factor 
D<== KATI0**2 
H<== 1-D 

Z<== PI*Scale_Faotor*Dbject»‘ 1*e-6 /Wavelength 
Do for I = 1 to 256 
Y< = Z*(I-1> 

If Y > 30 than 

RDATA (I) <== 0 

2 ]_ ss 

RDATA (I) <== (AIRY ( Y) -D*AIRY ( Y * R A T 1 0 ) ) *#2/H **2 
End If 
S<= I 
End Do 

(Record RDATA on separate fils) 

(Determines if a plot of the optics function is desired.) 
If OPTICS_FUNCTIDN_PLDT wanted then 
Ca 11 Plot (32, RDAT A, 1) 

End If 

(Calculates the line spread function of the optics 
function. ) 

Call LSF 

(Determines if a plot of the L3F of the optics function is 
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desir ed. ) 



If ?L3T_0F_LSF_0Pr ICS wanted then 
Call Plot (32, RDAr A, 1) 

End If 

Do for I = 257 to 512 



RD ATA (I) <== 3 
End Do 

(Dana are folded over again ana 

calculated .) 

Do for I = 1 to 256 

RD ATA (514-11 < == RDATA (I ) 



End Do 



: u r a a r transform 



RDATA (2 57) <== RDATA (256) 

(Record RDATA) 

(Fourier transform of the L3F of Optics is calculated.) 
Call FXFORH 

(Determines if plot of Fourier transform of Optics 
function is desired.) 

If PLOT_OF_FXFOR.i _ OF_D PTICS wanted then 
Call Plot (255 , RDATA, 2) 

End If 

(Calculation of the quotient of two Fourier transforms.) 
(This section takes the Fourier transform of the system 
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and divides it, point by point, by the Fourier transform 
of the optics. This yields taa modulation transfer 
function (MTF) of the a tmospa ar e. ) 

Input R DATA_OF_FXF ORM_OF_ SOURCE 
Do for I = 257 to 512 
RD ATA (I) < = = 0 
End Do 

Do for 1= 1 to 256 

RD ATA (514-1) C == R DATA (I) 

End Do 

Input RDATA_OF_FXF ORM_OF_OPTICS 
Do for I = 1 to 236 

RDATA (I) < = = I DATA (I) /RDATA (I) 

End Do 

(Determines if a plot of the dTF is desired.) 

If PL3 T_OF_MTF _OF _AT MO S_ wanted then 
Call Plot (32, RDATA, 1) 

End If 

(Calculation of C n 2 by curve fitting. 

This procedure is described in Chapter III.) 

(All "R" variables used in this calculation are 
initialized to 0.) 

R 17<==R 18< = = R1 9< = = R20<==R21< = = R22<==R2 3< = = R24<==R25<==0 
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(Analysis of data shows that computer tins is wasted 
beyond 96 points.) 

Do for I = 1 to 95 

(R18 - when summed is the (f(ii ) 10 / 3 . This is the 
total of the angular spatial frequencies squared.) 

R 1 9<== ( (I - 1 ► *Scale_Fao tor) **(10/3) 

R 1 7< = = R17 *■ R18 

(R19 - when summed is “.he (i(ii ) s / 3 . This is 
the total of the angular spatial frequencies.) 

R 1 9<== ((1-1) *Sca le_Factor) ** (5/3) 

R20<== HI 9 + R20 

(R21 - is the total number of points.) 

R2 1<== I 

Do While RDAT A (II > 3 

(R22 - when summed is the total cf the product of the 
natural logarithm of each point with R19.) 

R 22<== in(RDATA(I) ) * R 1 9 
R 23< = = R 22 + R23 

(R24 - when summed is the total of the natural logarithm 
of each po int . ) 

R 24< = = In (RDATA(I) ) 

R25<== R24 ♦ R25 
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End Do 



End Do 

(Outout matrix values foe debugging.) 

Output R17, R20, S21 ,R23, R25 
(Assign value to matrix. 

This section takes the above calculated values and 
sets up a matrix. The matrix equation solved is 
B = A- l *C.) 

(Cr eate 3 matrices A (2, 2) , 3 (2 , 1 ) ,C (2 , 1 ) . ) 



A (1, 


D 




R20 


A (1, 


2) 


< — 


R 21 


A (2, 


D 


< = 


R 1 7 


A (2, 


2) 


< “ 


R 20 


INVERT 


_ MAT RIX_A 


C (1, 


D 


< “ 


R 25 


C (2, 


1) 


< = 


R2 3 


B<== 


A * C 





(Output matrix values for debugging.) 

Output B (1 , 1) , B (1, 2) 

(Calculation of C n 2 using equation (2. 3) . ) 

R2 2<== E ( 1 , 1) /(-21 . 49* Range* R 20* (Wavelength** (-. 3333 3) ) ) 
Output »CNSQ=»,R22 

(The program now starts the prediction phase. It 
calculates a source assumed to have a Oaussian 
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distribution by meins of the computed source pattern 
below. Using similar Fourier transform techniques 
the program uses the calculated value- of C n 2 and predict 
the power incident on the target.) 

(Computed Source Pattern :A=AD*exp (-x**2/2*Sigma**2) 

(This is an "arbitrary' 1 Gaussian source pattern 
with a standard deviation for Sigma set equal to 2.) 

AO <= = A _Z E 3 0 

C<== S i gma*Scale_F actor 

Do for I = 1 to 25 6 

F<= (1-1) *Scale_Factor 

3 < = F**2/ (2 * C**2 ) 

Do While 3 > 13 

3 DAT A (I) <== 0 
End Do 

3D AT A (I) < = = AD *exp (-G) 

End Do 

(Record EDATA. This records the computed source.) 
(Determines if plot of computed source is desired.) 

If ?LOT_OF_ COMPUTE D_SD URCE_wanted then 
Call Plot (32, RDATA ,0) 

End If 

(Calculates the line spread function of tne computed 
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source. ) 



Cali LSF 

(Determines if a plot of the LSF of the computed source 
is desired.) 

if plot_of_lsf_of_comp_so(jr:e = I then 
Ca 11 Plot (32, RDATA ,0) 

End If 

(Caicuxarion or the Fourier crarsrorni or LSF or compared 
source. ) 

Do for I = 257 to 512 
R D ATA (I) < = = 0 
End Do 

Do for I = 1 to 256 

RD ATA (514-1) <== RDATA ( I ) 

End Do 

RDATA(2 57<== RDATA (256) 

Call FXFORM 

(Determines if a ploc of the Fourier transform of the 
LSF of the computed source is desired.) 

If PLOT _OF_FXFO RM _OF_EOMP(JT5 D_SO(JRCE_wa a t ed then 
Ca 11 Plot (255 ,RDATA,2) 

End If 

(Calculation cf the product of two Fourier transforms.) 
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(This is the Fourier zransfcriii of zhe computed source 
with the MTF cf the atmosphere.) 

Do for I = 1 to 256 
RD ATA (I) < = = 

RD AT A (I) _C03P UTED_S0U REE* RDATA (I) _FX FK_0F_I1TF 
En d Do 

(Product of FXF0RM5 of Source*Atmosphere*0?tics.) 

RDATA (!)<== 

RD AT A ( 1 ) _ PRODUCT_D F_2_FXFMS* RDATA (1) _FXFM_OF_OPTICS 
Do for I = 2 to 2 56 

RD ATA (I) <= = 

RD ATA (I) _PRDDUCT_OF_2_FXFKS*RDATA (I) _FXFX_OF_OPTICS 
RDATA (514-1) <==RDATA(I) 

End Do 

RDATA (2 57) <== RDAT A (256) 

(The above data predicts what the Fourier transform 

of the entire system is, using the calculated 

value of C n 2 for the atmospheric turbulence,) 

(Determines if a plot of the result of Fourier transform 

♦ 

is desired.) 

If PLOT OF FXFORfl PRODUCTS wanted then 



Ea 11 Plot (32, RDATA, 2) 



End If 
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(Inverse Fourier transform gives target LSF.) 

Call IN VERSE_FXF3R M 

(Determines if a plot of the inverse Fourier transform 
is desired.) 

If ?LOT_OF_INVERSE_FXFROM_w anted then 
Call Plot ( 32 1 R D AT A , 0 > 

End If 

(Converts target's one dimensional LS? to a 
two dimensional PSF by Abel transform.) 

Call ABEL 

(Determines if a plot of the Aoel transform is desired.) 

If P LOT _OF_ ABEL _X F ORM_ wan ted then 
Call Plot (32 f RDAIA, 0) 

End If 

(Calculates the fraction of power inside circle of radius 
R. This predicts the fraction of power that will be 
incident on target. ) 

RD AT A ( 1 ) <= = 0. 25*PI*RDATA(1) 

Do for I = 2 to 256 

R D AT A (I) < = = 2*PI*RDATA (I) +RDATA (1-1) 

End Do 

(Determines if a plot of the fraction of power is desired.) 
If PLOT OF POWER wanted then 
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Call Plot (32, RDAT A, 0| 

End If 
End Laser 

(The following subroutines ire used in the Algorithm 
Laser . ) 

S ubr out ine_Tr ansf ar ( F ile_ Number I 

(The transfer subroutine aets raw data, booh Dreaatble 
and data from the waveform massage sent by the 
465. It processes t.ha message by finding the minimum 
value of the array then subtracting- this value from 
each element in the array. This ’'zeros" the array.) 

Output 'Ensure Equipment set up properly' 

Output 'Continue when ready' 

S<== (-1) 

Do while S = -1 

("DATA" is the I/O buffer where the data from the 468 are 
sent. The status of the buffer is read while daza are being 
transferred. When the transfar is compleza, the interface 
is cleared.) 

Data<== Trans ferrea_Data 
If Data_Trans fer_Complate then 
E <== 0 
End If 
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Cl ear_ Inter face 
N<== 1 

(This is where healer format stops and data begin.) 

A<== po s(Daf a, * %* ) 

Do for I = A+2 to 687 by 16 
Do for J = 1 to 16 

If J + I>5 88 t hen 
RJ<== 0 
N< = = N+ 1 ' 

End If 

RJ<== na m (Da ta (It J) ) 

(The numerical values of each element are stored in YD A TA. ) 
Y DAT A (M| <== RJ 
N<== Htl 
End Do 

End Do 

E<== minimum (YD AT A) 

Do for I = 1 to 51 2 

YDATA (I) <== YDA I A (I) - E 
RD AT A (I) <== Y DAT A (I| 

End Do 

Store_3 DAI A (Fil e_M umber ) 
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Return 



End Transfer 

(This subroutine calculates the Line Spread Function 
for the previously recorded point spread function 
by using equation (1.1). > 

Subrout ine_LSF 

(Depending on which file number has been passed down 
from the calling suocoutine, either scale or source 
data are loaded.) 

If File_Number = 3cale_Data then 
Input Scale_Data 
RD ATA<==SCALE_DATA 

Else 

Input RDAIA 
End If 

Do for I = 1 to 512 
ID ATA (I) <== 3 
End Do 

Do for I = 1 to 256 

ID ATA (I) 0= S DATA (I) 

End Do 

(Plots have shown that computer time is lost and no 
valuable information is gained beyond about 24 points.) 
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Do for I 



1 to 24 



J<== 1 
Output I 
Q<= IDATA (I| 

Do While R < 24 

R<== SORT (I*I+J* J| 

Q<== 

Q + 2* ( (1 - fraction (3 ) ) *1 (int (3) ) + fraction (?) *1 (in t ( R) + 1) ) 
J<== J> 1 
End Do 

IDATA(I)< = = 0 
End Do 

Do for I = 1 to 2 % 

RD ATA (I) < = = I DATA (I) 

End Do 

If Fiie_Number = Scaie_Data rhen 
Store_Scaie_D ata 

Else 

S t ore_RDAT A 
End If 
Return 
End LSF 
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(This subroutine calculates the Fourier transform 
of the given data using the Co o ley-Tukey Algorithm. 
If the inverse statement is true, then the inverse 
Fourier transform is calculated.) 

S ubrout ine_FXFORM 
Set_Rad ian_Mode 

(2**9=512 which is the number of points.) 

N<== 9 

Do for I = 1 to 51 2 
ID ATA (I) < = = 3 
End Do 

If File_Number = Scale_Data then 
Input Scale_3ata 

Else 

Input RDAIA 
End If 

T<== FI/2** (N-1) 

Do for J = 3 to 2 ** (M- 1 ) - 1 
(31 is a bit inversion subroutine.) 

Call 31 (J,P, N-1) 
c<== cos (P*r) 

If INVER SE.FXFORM = True then 
Fig 7<== 1 
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Fig 7<== 0 



End If 

P<== sin (P*r)*(1-2*Flg7) 

Do for I = 2*R0*J+1 to 2*R0* J + RO 
R1<== RDATA ( I ) 

R2<== RDATA (I + RO) 

33 < == IDATA (I) 

R4<= = IDATA (I + RO ) 

RD A TA (I) <== R1+R2*C+R4*P 
IDA TA (I) < = = R 3 + R4 *C-R2*P 
RDA TA (I + RO) < = = R1-R2-C-R4*? 
IDATA (I + RO) < = = R3-R4*C+R2*P 
End Do 
End Do 
Output M 

End Do 

(This section is for re-ordering the block.) 

Do for I = 0 to 2**N-1 
Call 3I(I r J, S) 

Do While I-J < or = 3 
Do While I * J 

?<== RDATA (1+ 1) /SQRT (2**N) 
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Z<= = IDATA (1+1) /SQRT(2**N) 

RDATA (1 + 1) < = = RDATA (J + 1 ) 

IDA TA (J + 1) < = = IDATA (J + 1 ) 

RDATA (J + 1) < = = P 
IDATA (J+1) < = = Z 
End Do 

RDATA (1+ 1) <== RDATA (1+1) /SQRI (2**N) 

IDATA (I + 1) <== IDATA (1 + 1) /S QR T ( 2 ** N) 

E nd Do 

If File_Nu:nb er = Scala_Da-a then 
Stcre_Scale_Data (File_Number) 

E lse 

3tors_RD AT A ( File_Siu m ber) 

End If 

End Do 

Set_Deg rees _ Mode 
Return 
End FXFORM 

(This subroutine takes the binary number PI containing 
P3 bits and inverts it end for and, e.g. 310111 becomes 
111010.) 

Subroutine BI(P1,?2,?3) 

? 2 <== 0 
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pu <== pi 



Do for Z = 1 t: ?3 
P 4 <== P4/2 
P 2 <== 2*P2 

If fraction(P4) ir 0 than 
? 2<== ? 2 + 1 
End If 

P 4 <== in* ( P4 1 
End Do 
Return 
End 31 

(This function calculates Airy ( x) =2*J 1 (xj /x where J1(x) 
is the Bessel function cf order one.) 

Function AIRY(PI) 

If PI < 0 then 

Output 'ERROR -AR3 UMENI LESS THAN O' 

Step 
End If 

If PI = 0 then 
R 4 <== 1 

Return R4 
End If 
R5<== 0 
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If PI > 15 then 



R 6 <== 90 f PI/2 
3 1 2< = = 1 . 4 *P 1 + 60/PI 

End If 

If PI < 5 then 

R 6 <= = 20 *■ 10 *P1-P1**2/3 
R 1 2< = = 6 ♦ PI 

Else 

R 6 <== 20 <■ 10*P1-P1**2/3 
R 1 2<== 1. 4*P1 + 5 0/P1 

End If 

R1 2<== maxi mum (in t (P 12) , int (3 + P 1/4) ) 
Do for LI = R12 no R6 by 3 
R 8 <== 1*9-23 
R 1 3<== R 14<== 0 
If M/2 = inn ( M/2) then 
Flgl Q<== 0 

Else 

Flgl 0 < = = 1 
End If 

Do for J = 1 to 3-2 

R 1 5< = = 2 * (M-J) *R3/P1 -R13 
R 1 3< = = R 8 
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38 <= = 31 5 



If M-J-2 = 0 than 

R4 < == a is 

End If 

If FlglO = 3 than 
Fig 10<== 1 

Else 

Fig 1 0<= = 0 
End If 

R 1 4< = = a 14+2*R8*FLg 10 
End Do 

a 1 5< = = 2*R 8/? 1-R1 3 
HI 4< = = E14+R1 5 
R4 <== R4/R1 4 

If (abs (R4-R5) -abs(R'4*i *e-6) ) < or = 

R4<= = 2 * R4/P 1 
End If 
Return R4 
E 5 <= = R4 

» 

End Do 

Output 'ACCURACY SOT DETAINED' 

Return R4 
End AIRY 



0 “hen 
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(The subroutine Abel takes a 01 e-dimensionai 
line spread functiou and calculates a two 
dimensional point spread function.) 

Subroutine ABEL (File_Number) 

N<== RD ATA (1) 

RD ATA ( 1)<== 1. 4 *R D ATA ( 1)-1.3*RDATA(2)+. 4 * ED ATA (3) 
Do for I = 2 to 54 
H<= R D ATA (I| 

RD ATA (I) < = = . 4*N+.2*M-. 5*RDATA (1+1) 

N < = M 
End Do 

Do for I = 1 to 54 

RD ATA (I) < = = 3 DATA (I) /(2*SQRT ( (I + . 1) **2-1*1) ) 
Do for J = 1 to 54 
R DATA (I) <= = 

RDATA (I) +RD ATA ( J) /S QR T ( ( J+ . 1 ) **2-1*1) 

End Do 

RDATA (I)<= = RDATA (I) /PI 
Output I 
End Do 
Return 
End ABEL 
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(This subroutine prints out data in tabular fori.) 

S ubrout in e_T ABLE (Fils _Nuaber) 

Input R DATA (File_>I umber) 

Do for I = 1 to 32 

Do for J = 1 to 15 
Output RDATA ( 16 (1-1) *J| 

End Do 

Output RDATA (16*1) 

Return 
End TABLE 

(This subroutine plots the desired data.) 

(Each time a plot is called, tiis subroutine 
plots the particular graph and labels it. 

The horizontal axis is labeled according to the first 
and last parameters passed by the calling subroutine.) 
Subroutine_PLOT (Pi ,Fi le_N umber , ?3) 

Input RDATA (File_N umber) 

A<== 0 
B<== PI 

C < — — minimum (RDATA ) 

D<== maximum (RDATA) 

XM IN <== -. 1 * (B-A) 
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XMAX<== B+.05*(3-A) 



YMIN<== C-. 1* (D-C> 

YMAX<== D+. 05* (D-T) 

E<== B 
F<== 10 

(Using the value of PI, the horizontal increment fc 
plotting and labeling is determined. ) 

If ?1 = 512 then 
3 < = 64 

Else If PI = 256 then 
3<== 3 2 

Else If PI = 64 then 
3 <= 8 

Else 

3 < == 4 
End If 
Plot B , C, 1 

(This lines and places tic marks on the horizontal axis.) 
Do for I = E to 0 by - G 

Plot I* (B-A) /E+A, C, 2 
Plot I* (B-A) /E+A, C+ (D-E) / 150 , 2 
Plot I* (B-A) / E+ A, C , 2 
End Do 
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Do for I = 0 to ? 



Plot A,I*(D-::)/F+C,2 
Plot A + (B- A ) /I 5 0 , 1 * (D- C) /F+- T , 2 
Plot A , I* (D-C) /Ft-Z ,2 
End Do 

Charact er_Size 1.2,1,. 7,0 

(This numerically labels the /erzical axis.) 

Do for I = F to 0 by -1 

Plot A-.075* (B-A) ,1* (D-T) /F + C,1 
Label I/F 
End Do 

(This sets up a character string for labeling the vertical 
axis. ) 

If P3 = 0 then 

AS TRI NG<== 'DATA POINTS" 

L<= 1 

Elsa If P3 =2 then 

A S TRING<== • LINES/MI3RD RADIAN* 

L< — 1 / (2*L) 

End If 

(This numerically labels the horizontal axis.) 

Do for I = 0 to E by G 
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Plot A + (I/E-. 025) * ( B- A* , C - . 0 25* (D-C) , 1 
Label I*L 
End Do 

If P3 # 0 then 

L<= Scale_Factor 
End If 

(This labels the horizontal axis.) 

Plot . 4 *(B-A) + A 0 5* (D-C) +C, 1 
Label A STRING 

Plot 07* (B-A) + A, .3*(D-C) + C, 1 
Charact er_S ize = 1.2,1, .7, 90 
ASTR IN G <== ’NORMALIZED INTENSITY' 

Label A STRING 

Charact er_Size = .5, 1,1. 5,0 
(This plots data.) 

I<== 0 

Plot I , maximum (RD A TA (I ) ) , 1 
Do for I = 1 to PI 

Plot 1-1, S DATA (I) 

End Do 

If P 3 = 0 then 

Input PLOT_LABEL 
IS THING <== PLOT LABEL 
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End If 



Plot . 6*B, . 9*D, 1 
Label • ELOT OF' 
Plot . 6*B # . 87*D, 
Label ISTRING 
Plot . 6 *B, . 84*D , 
Label ESTRING, • 
Return 
End PLOT 



1 

1 

Laser ' 
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V. 



CONCLUSION 



The work reported in this thesis supports the model 
predicted and measured in Crittenden, and ethers, for the 
long exposure case [Ref. 2]. The vidicon, in replacing the 
mechanical slit scanning system, shows no degrading of the 
signal data. It also demonstrates that a good approximation 
for the point spread function may be made by recording a 
single TV line through a laser spot. This line is then used 
to calculate the one-dimensional line spread function. The 
linearity of the video tape recorder is seen in the results 
of the MTF and C n 2 . The 463 oscilloscope is the workhorse 
for the entire system. It effectively displays, stores, 
digitizes, and transfers data on a real time basis. 

The overall system does not have the capability of the 
system described by Crager. However, the relative 
simplicity cf the structured programming technique coupled 
with a lower equipment cost demonstrates that comparable 
measurements and data evaluation using this equipment can be 
made. Further, investigation should include use of the Data 
Precision 6000 digital waveform analyzer thereby allowing 
data to be sampled and analyzed on a near real time 
basis. 
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APPENDIX A 



0: "HE-NE : SOURCE ON 9 , SCALE ON /,& COLLIMATED ON 5": 

1: "GA-AS: SOURCE ON 8 , SCALE ON 6 ,&CCLLIiMATED ON 4" : 

2: ent "SELECT CODE FOR PRINTER" ,A;aev "print", A 

3: tint 1 , z , c ;tmt :2 , t8 .1 ,z; tint 3,t8.1 

4: a im I$.[1024] ,R[512] ,Y[ 512] ,C$(2G] ,D$[20] ,E$[20] 

5: aim I [ 512 ] ,A$ [ 32] , B$ [2 0] ,A [ 2 , 2 ] ,B ( 2 , 1] , C ( 2 , 1] 

6: but "LATA" ,1$ ,.3; ina l ,1 ,R;0*F*W-*-V+U; 1*L; but "I",I$,3 

7: beep; asp "EN1EE LASER USED";wait 15U0 

8: ent "HE-NE=1 ,GA-AS=G" ,W 

9: it v\=l; "riE-NE " -*-E$ ; 9 +C; 3it'P 2 

10: "GA-AS"*E$; Q+Q 

11: beep; asp E$&" SOURCE DATA RECORDED?" ;vvait 1500 
12: ent "1=YES , 0=NO" ,V 

13: it V = 0 and to=l; stg 0;cll 'TRANSFER ' (0) ; j rap 2 

14: if V=0 and to=G;srg 0 ;cll "TRANSFER ' (Q) 

15: beep;ent "PLOT OF SOURCE? , 1= YES , 0 = NC" ,R 
16: it R=l;Cll 'PLOT ( 512 ,Q , G) 

17: beep; asp E$&" SCALE DATA RECORDED? "; wa it 1500 
18: beep; ent " 1=YES , G=NC" , U 

19: ir U=0 ana W=l;sfg 0 ; ell ' 'TRANSFER '( 7-*Q) ; jmp 2 
20: if U=C ana W = 0; sfg G;cll " 'TRANSFER '( 6+0) 

21: beep ;asp "WHAT IS SCALE FACTOR?" ; uai t 1000 

22: ent "SCALE FACTOR=?;G GETS PLOT" ,S 

23: it S = C ana U i 0 ; v 2 -*• v ; c 1 1 ' P LOT '(512,0,0) 

24: if S=0 ana U = 0 ; ell ' 'PLOT ' ( 5 12 ,Q , 0) 

2 5: it S=0; jmp .-3 

26: S ♦L 

27: it W = l; 9-*-.T-*-0;j^mp 2 

28 : 8-.T^0 

29: lat T,I$ ,R[*] 

.30: pos ( 1$ , "R: " ) -B 
.31: val(I$[B+2,B+4] )-rl 
.32: pos (I "IT:".) -6 
33: I$.[B+.3,E+4]-A$ 

34: lot .T'-2 , 1 $ ,R [ * ] 

35: pcs ( I $ , " R : " ) ♦B 
.36: val(IS[B+2,B+4].)*r2 
37: pos(I$,"LT:"0-B 
38: I $. [B + 3 , B+4 ] -*-B$ 
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39: let .T-2,I$,R[*] 

40: pos ( 1$ , "R: " ) -E 
41: val(I$[B+2,B+4] ) -*• r 3 
42: pos(I$,"LT:'\)-B 
43: I $ { B+ 3 , E+4 ] -m3 $ 

44: if r2 = rl and AS=B$ ana r2=r3 ana B$.=C$;jmp 2 

45: beep;asp "ERROR IN TIKE SCALE" ; stp 

46: "B" : ldf Q, I $ ,R [* ] ; ina Y;ara R + Y;C+A+B-*C 

47: for 1=1 tc 512 

48 : A+Y [ I ] ♦A ; if I < 51 ;B+Y [ I ] ♦B 

49: if I >4 63 ; B+ Y [ I ] ♦B 

50: next I;prt A,B 

51: B/100-D; A-512*D^A 

52: tor 1 = 1 tc 512; Y [ I ] -D- Y [ I j ; C+Y [ I J *C 
53: it C < A/2; I-J 
54: next I;prt C,J 
55: tor K=1 tc 256 

56: Y [ ( J.-K+L) mca512 + l] + Y [ (K+J-l) ir.cc512 + l] -*■! [ R ] 

57: next K;ina R;ara I-»R 

58: "IMAGE POINT SPREAD FCN"-I$;rct JG,I$,R[*] 

59: beep;ent "PSF CU.TPUT = 1 AND/OR CONI" ,Z 
60: beep ; ent "PSF PLGT = 1 AND/OR CONT" ,Y 
61: if Z=l;cll "TABLE" (10) 

62: if Y = l;cli "PLOT" (32 ,10,1) 

63: min (I [ * ] )*r22;raax(I [*].) *r23 
64: prt r22,r23 

65: Deep; ent "LSF CUTPUT. = i AND/CR CCNT" ,Z 
66: beep; ent "LSF PLOT=l AND/OR CONT" ,Y 

67: ell "LSF"; "IMAGE LINE SPREAD FCN " -I $ ; ref 1G,I$,R[*J 
68: it Z = l;cll ".TABLE" (10) 

69: if Y = l;cll "FLCT"(32 ,10,1) 

70: "CALCULATION OF FXFORM OF LSF":laf 10 , 1 $ , R [ * ] ; ina I 
71: ara R-I 

72: I. [ 1] ♦R [1] ; f or 1 = 2 tc 256 
73: I[I]+R[I.]-*-R[514-I] ;next I 
74: R [ 25 6 ] +R [ 257 ] ;rct 10,I$,R[*] 

75: ell "FXFORM'; "FXFORM OF IMAGE LSF"-I$;rcf 10,I$,R[*J 
76: Deep ; ent "FXFORM OF LSF OUTPUT. = l AND/CR CONT" ,Z 
77: beep;ent "FXFORM OF LSF PLOT = l AND/CR CCNT" ,Y 
78: if Z = 1 ; ell ".TABLE" (10) 

79: if Y=l;cll "PLOT" (32 ,10,2) 

80: " B=D I A OF CBSCUR/DIA OF OBJ LENS IN METERS": 

81: "B=0 . 064/ . 16 4 " : 

82: " r 1=SCALE OF DATA IN MICRORADIANS , rl-1 . 60" : 

83: "0=D IA OF OBJ LENS IN METERS ;O=0 . 1 64M" : 

84: " W =WA V ELE NGT'H IN METERS": 
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85: "CALC CF DIFFRACTION LI NIT POINT SPREAD FCN" : 

bo: ina R;if W=l; 6 . 328e-7-»-W ; jmp 2 
8 7: 9.0 5e- 7-*-W 

88: . 064/. 164-E; . 164+0 ;L^rl 

8S: 2*D; 1-D-H; n*rl*0*le- 6/w-Z 

90: for 1=1 to 256 

91: Z* ( I -1 ) ♦Y; if Y>30;gtc +3 

92: ("AIRY ' (Y) -D* 'AIRY ' (Y*B) ) ~2/H~2-R[I] 

93: I-*-K;gto +2 
94: 0-R[I] 

95: txa .5; asp R[I];next I;"0PTICS FUNCTION "♦I $ 

96: ref 11,I$,R[*] 

97: beep; ent "OPTICS FCN PLOT = l AND/OR CONT" ,Y 
98: if Y = l;cll "PLOT ' (32 ,11,1) 

99: sfg 5 ; c 11 ’'LSF ' ; " LSF OF OPTICS FCN " -I $ ; ref 11 , 1 $ , R [ * ] 

IOC: beep ;ent "PLOT OF LSF OF OPTICS FCN = 1 AND/OR CONI" ,Y 
101: ir Y=l;cll "PLOT '(32,11, 1) 

102: sfg 5?ina I;ara R-*I; I. [ 1] ->-R[ 1] 

103: for 1=2 to 256; I [ I ] *R [I ] -R[ 514-1] ; next I 
104: R [25b ] ♦R [ 25 7 ] ; r cr ]1 ,1$ ,R[*j 

105: ell "FXFORM V'FXFCRM OF LSF OF CPT'I CS" -*-1 8 ; r.cr 11,I$,R[*] 
106: beep ;ent "PLOT OF FXFM OF OPTICS=l AND/OR CONT" ,Y 
107: if Y = l;cll * 'PLOT ' ( 256 , 11, 2 ) 

108: "CALCULATION CF QUOTIENT CF TWO FOURIER TRANSFORMS": 

109: lar 10 ,I$,R(*] ; ina I;ara R+I 
110: idf 11 , 1 $ , R [ * ] 

111: for '1=1 to 256 

112: I [I ]/R(I] -R[I] ;next I 

113: "M1F CF SYSTEM" -I $ ; re f 10,1 ?,R(*] 

114: beep ; ent "PLOT OF MTF OF SYS=1 AND/OR CONT" ,Y 
115: it Y=l;cll "PLOT ' (32, 10 , i) 

116: 0->rl7-^rl8^rl9-*-r20-»-r21-»r22-*-r2 3>r24^r2 5 

117: tor 1=1 to 96 

118: ( (I-1)*L) *(10/3) -rl8 

119: rl8+rl7-*r 17 

120: ( (I-1)*L) * (5/3) -rl9 

121: rl S+r20-*-r 20 

122: I-*-r21 

123: it R(I]<=0;gtc "CC" 

124: (In ( R [ I ] ) +r24) *rl9+r23-r23 

125: r24+r 25-*-r25 

126: "CC" : next I 

127: prt r 1 7 , r 20 , r2*l , r 2 3 , r2 5 

128: r 20-*-A [ 1 , 1] ;r 21->A[1 ,2] ; rl7*A[2 4 1] ; r2 0+A [ 2 , 2 ] 

129: r25-*-C [ 1 , 1] ; r23-*C [ 2 , 1] 

130: inv A-*-A;mat A*C-*B;tlt 5;aprt B 
131: "CALCULATION OF CNSQ" : 

132: B [ 1 ,1 ] / (-21.4S*14 5*r20*W*(.-. 3 3 33 3) ) -r2 2 
133: prt " CNSQ=" ,r22; stp 
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134: cfg ; laf 2 

135: "TRANSFER" : beep; asp "ENSURE EQUIP. SET UP PROPERLY" 

136: wait 1500 

137: asp "PRESS CONTINUE WHEN READY " ; stp 
138: dsp '"wait ing " ; wai t 1500 
139: buf "DATA" ; 0-Z 
140: tfr 703," DATA" 

141: asp "transferring" 

142: ras ( "DATA" ) -E; it E=-l;jrcp 0 
143: clr 703 

144: dsp "setting output" 

145: pos (1$ , "%") -A; 1-N 

146: for I=A+2 tc 687 by 16 

147: for J=1 tc I6;if J + 1 >68 8 ; C-»-r J ; jmp 3 

148: num(I$[I+J] ) *rJ 

149: r J -*• Y [ N ] 

150 : N + l-*-N 
151: next J 
152: next I 

153: min ( Y (*] ) +E; for 1 = 1 cc 512; Y [ I ] -E-*Y [ I ] ; next I 
154: ina R;ara Y->-R; rcr pi , I $ ,R [ * ] ; r et 
155: "LSF":if tig 5; laf 11 ,1 $ ,R[ *] ; jmp 2 

156: laf 10,I$,R[*] 

157: ina I;ara R-*-I;tcr L=1 tc 24 
158: 1-J;csp I ; I [ I ] -C 
159: y ( I * I+J* J ) -*-R 

160: 2* ( (1-rrc (R) ) * I [ mt(R) ]+frc(F.)*I[ int (R) +1] ) +q-*-q 
161: J + l-*-J;if R<24;jmp -2 
162: Q-L[I] 

163: next I;ina R;ara I-»R; if rlg5;rcf 11 , 1 $ , R [ *.] ;crg ;ret pi 
164: ret pi 

16 5: "FXFORM" :rad;S-*-N; ina I; if rlg5;lcf 11 , 1 $ , R [ * ] ; jmp 2 

166: laf 10,I$,R{*] 

167: tt/ 2~ (N-l) *T 

168: for M=1 tc N ; 2" (N.-Pi) -rO 

169: for J = 0 tc 2 * (M-l).-l; ell *'BI ' ( J , P ,N-1 ) 

17 0: cos (P*T) -K3 ;sin (P*T) * ( 1-2 *tlg 7) -»-P 
171: for 1 = 2* r 0*J+1 tc 2*r0*J+rG 
172: R[I]-rl;R(I + rO]-*r2 

173: I [ I ] ♦r 3 ; I [ I + r 0 ] -*■ r 4 

174: rl+r 2*C+r 4*P-*-R ( I ] ; r 3+r4*C -r2*P*I [ I] 

17 5: rl-r2*C-r4*P-R[ 1 + rC ] ; r 3-r 4*C + r 2* P- I [ I + r0 ] 

176: next I; next J;asp M; next M 
177: for 1 = 0 to 2~N-1 ; ell ' '61 ' ( I , J ,N ) 

178: if I-J >0;gto "£B" 

179: if I=J ;g to "INC" 

180: R[I+1]/V(2"N)^P;I[I+1]/V (2*N)-»Z 
181: R[J+1]*R[I+1] ;I [ J + l] *I[I+1] 

182: P-R[J+1] ;Z-I [J+l] 
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163: "INC" :R{I + 1] /V (2~i\') -R(I+1] ; I [1 + 1] /v ( 2~M + 1 [ 1 + 1] 
164 : "EB" : nex c I 

185: u rlg5;rcr 11 , 1 6 ,R [ * ] ; a eg ; re c 

166: aeg;rec 

187: "BI" :0-»p2;pl-p4 

188: lor Z=1 tc p3 

18S : p4/2-*-p4 ;2*p2-*-p2 

190: it r rc ( p4) #0; p2+l-*-p2 

191 : int ( p 4 ) -*-p4 

192: next Z ; ret 

193: "AIRY" : if pl< 0 ; beep ;asp "er rcr-argunient<0" ; s tp 

194: it pl=0; Or 4; re t r4 

195: iOr5;if pl>15;jmp 2 

196: 2G+ 10*pl-pl”' 2/3-*-r6 ; jrcp .2 

197: 9U + pl/ 2-*-r6 

198: ir pl< 5; 6 + pOr 12; jrrp 2 

195: 1 . 4*pl+6 C/pOr 12 

200: rr.ax(int irl2) ,.int(3+pl/4) Or 12 

201: tor R=rl2 tc r6 by 3 ;le-2o -»r6; C-*-rl3-»r 14 

202: srg 10; it ¥/ 2=int (M/ 2) ;crg lo 

203: ter J=1 tc K-2 ; 2* (JOJ ) * r 8/ pi- rl 3-*-r 1 5 ; rb-*-r 13 

204: rl5-*-r8;ir M-J-2 = G; r 15-*r4 

205: cir.f 10; r 14 +2* r 8*1 lg 1 0-*-r 14 ; nex t J 

206 : 2*rfcs/pl-rl3->rl5 

207: rl4+r 15-»r 14; r4/rl4-*-r4 

206 : ir abs (r 4- r5 ).-so£ (r 4* le-6) <=G; 2* r4/pl--r4 ; r et r4 
209 : r4-*-r5;next to 

210: ceep;dsp "ACCURACY NOT OBTAINED'" ; wait 1500;ret r 4 

211: "1AELE " : idl pi,I$,R[*J 

212: tor 1=1 tc 32;icr J=1 tc 15 

213: wrt "print. 2" ,R [ 16 (I -1) +0 ] ; next J 

214: wrt "pr int . 3 " ,R [ lb I ] ; nex t I; ret 02 

215: "PLOT": lot p2,I$,R[*J 

216 : 0-*-A;pO5;m.in(R( *] ) -*-C ; rr.a x (R [ * ] ) -<-D 

217: SCI A-.l(B-A) ,B-h.C5 «£-A) ,C-.1(D-C) ,0-h. 0 5 (D-C) 

218: JOE ; 1C+F 

219: it pl= 51 2; 6 4-*-G 

220: if pl=256 ; 32+G 

221 : 1 1 pl = 6 4 ; 6+G 

222: it pl=32;4^G 

223: pit B , C / 1 

224: ror I=E tc 0 by -G 

225: pit I (E-A)/E+A ,C ,2 

226: pit I ( E-A) /E+A ,C+ (D.-C) / 150 , 2 

227: pit I (B-A) /E+A ,C,2 

228: next I 

229: tor 1=0 tc F 

230 : pit A, I ( D-C) /F+C , 2 

231: pit A+(B-A) /150 , I (D-C) /F+C ,2 



93 



222: pit A, I (C-C) /F+C, 2 
253: next I; pen 
2 34 : c siz 1. 2 , 1 ,. . 7 
235: fxa 1 

236: for I=F tc 0 by -1 

257: pit A- . 07 5 (B-A) , I (D-C)/F+C , 1 

258: lbl I/F 

25S: next I 

240: fxa 0 

241: it p3=G; "DATA POINTS" -A$ ; 1-L 

242: if p3=2; "LINES/MICRCRADIAN" ♦A$ ; 1/ ( 2*L) -*L 

243: for 1=0 tc E by G 

244: pit A + ( I/E- . 0 25 l) (B-A) ,C-.G25(D-C) ,1 
245: lbl I* L 

246: next I; if p3*G ana p3*2; ,, MICRCRADIAKS ,, -A$ 
24 7 : if p3 = 2; S*L 

248: pit . 4 ( E-A.) +A ,>-. 0 5 ( D-C ) +C ; lbl A$ 

24 S : pit 07 (B-A) +A , ..3(D-C)+C; csiz 1.2,1, .7,50 

250: "NORMAL I ZEE INTENSITY"-*^.-?; lbl A$ 

251: csiz .5,l,1.5,G;0-*-I 
252: pit I ,rnax(E[*) ) ,1 

253: tor 1=1 tc pl;plt I -1 ,R [ I ] ; nex t I 
254 : csiz 1. 2 , 1 . 7 , 0 

25 5: if p5#0,*gtc +2 

256: beep ;ent "PLOT LABEL?" ,1$ 

257: pit . 6B, . SD , 1; lbl "PLOT OF" 

258: pit . 6B , . 8 7D , 1 ; lbl 1$ 

259: pit . bS , . 84D , 1; lbl ES," LASEF" 

26C: pen,-cfg ;ret 0-*Y 
*28625 



94 



0: "COMPLIED SOURCE PATTERN": 

1: ina R; 1000-^A ;,2*LC ; tr k l;faf 10 

2: ror 1=1 tc 256 

3: L*(L-J)-F 

4: F ~2/ ( 2*C~ 2 ) -*-G 

5: it G>13;jmp 2 

6 : A*exp(.-G) +R[ I] ;I-K; jrcp 2 

7: 0 -*-R [ I ] 

8: next I 

S: "COMPUTED SOURCE" -I $ 

10: ret 10,IS,R[*] 

11: beep ;ent "PLOT OF COMPLIED SOURCE = l AND/OR CONT" ,Y 
12: it Y = l;cll "PLOT" (32 ,10, 0) 

13: " CONVERTS PSF TO LSF " : 

14: cl 1 'LSF' 

15: "LSF OF COMPUTED SOURCE" -1 5 

16: ref 10 , 1 $ , R [ * ] 

17: Deep ; ent "PLOT OF LSF OF COMP SOURCE = l AND/OR CONI " ,Y 
18: if Y = l;cll "PLOT" (32,10,0) 

IS: "CALCULATION OF FXFORM OF LSF OF COMPUTED SOURCE": 

20: ina I;ara R-*-I 
21: I [ 1 ] -*• R [ 1] 

22: tor 1=2 tc 256 

23: L[I]+R[I]+R[514-I] ;next I 

24: R [ 256 ] ->-R [ 257 ] ; ret 10,I$,A[*] 

25: ell' 'FXFORM' 

26: "FXFORM OF LSF OF CALC SOURCE " ■*! 8 ; ret 1G,I$,R(*] 

27: Deep ;ent "PLOT OF FXFM OF CALC SRCE = 1 AND/OR CCNT" ,Y 
28: if Y = l;cll "'PLOT '(256,10,2) 

2S : "CALCULATION OF PRODUCT OF 2 FX FORMS " : 

30: laf 10 , 1 $ , R ( * ] ;ina I;ara R+I;trK 0 
31: let 1 0 , 1 $ , R [ * ] 

32: for 1=1 tc 256 

33: R ( I ] *1[I]-R[I] ;next I 

34: ref 1G,I$,R[*] 

35: "TRANSFER FCN OF SOURCE* ATMOSPHERE*ORTICS " : 

36: lat 11,1$ ,R(*.J ; ina I;ara R-*I 

37: III] *R[1]*R[1] 

38: for 1=2 tc 256 
3 S: R ( I ] * 1 [ I ] -*-R [ I ] -*-R [514-1] 

40: next I 

41: R[256] ♦R [ 257] 

42: "FXFM OF SGURCE*ATMOS*OPTICS" -I $ 

43: ref 10,I$,R[*] 

44: beep;ent "PLOT OF FXFM PRODUCTS = l AND/OR CONT" ,Y 
45: it Y = l;cll "PLOT" (32 ,10,2) 
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46: "INVERSE FXFORM GIVES TARGET LSF " : 

47: sf g 7 ; c 1 i ‘ 'F XFCRM ' ; cf g 7;"INVERSE FXFCRM"-* 1$ 

48: ref 10 , 1 $ , R [ * ] 

4S: oeep;ent "PLOT OF INV FXFM-1 AND/OR CCNT" , Y 

50: if Y=1 ; c 11 ' 'PLOT ' ( 32 ,10 , 0) 

51: "CONVERTS LSF TO PSF BY AEEL TRANSFORM": 

52: Cll "ABEL "ABEL XFOP.M PSF"-I$ 

53: ret 10 , 1$ ,R[*] ; sfg 0 

54: beep ;ent "PLOT OF ABEL XFORM=l AND/OR COKT",Y 
55: if Y=1 yell ‘ 'PLOT ' ( 32 , 10 , 0) ;stp 

56: "CALCULATES FRACTION OF POWER INSIDE CIRCLE OF RADIUS R" 

57: .25*u*R[l]*R[l] 

58: for 1=2 tc 256 

5S : 2tt* I * R ( I ] +R [ L-l] -*-R ( I ] 

60: next I 

61: "POWER FRCTN INSIDE CIRCLE" >1$ 

62: ref 10,I$,R[*] 

63: ceepjent "PLOT OF POWER = l AND/CR CCNT" ,Y 
64: it Y=l;cil ' 'PLOT ' ( 32 ,10 , G) 

65: stp 

66: "LSF": it tlg5;lof 11 ,1 $ ,R [ * J ;nir.p 2 

6 7 : lof 10 , 1 $ , R [ * ] 

68: ina I;ara R-*-I;tor 1 = 1 tc 24 
6 3: 1-J;cisp I ; I [ I ] -Q 

70: y ( I * I+J* J ) ♦R 

71: 2* ( (1-rrc (R) ) *l [ int ( R) ] +trc (R)*I. [mt <R) +1] ) 

72: J + i-*-J;if R<24;gir,p -2 
73: 0-1(1] 

74: next I; ina R;ars I-R; if rlgSjrcr 1 1 , 1 $ , R [ * ] ; c tg ;ret cl 
75: ret pi 

76: "FXFORM" : raa; S-N; ina I;it flg5;lar 11 , 1 $ ,R [ * ] ; imp 2 

77: let 10,I$,R[*] 

78: tt/2~(N-1)VT 

73: fer M = 1 tc N; 2“ (N-M) ♦rO 

80: for J = 0 tc (K-l).-l ; ell ' '31 ' ( J ,P ,N-1) 

81: .ccs (P*T) -C ;sin ( P*T) * (1-2* tig 7 ) -P 
82: fer I=2*rO*J+l tc 2*rG*J+rO 
83: R[l]-rl;R[I+r0] -r2 
84: I [ I ] -r 3; I. [ I + r 0 ] -r 4 

85: rl+r2*C+r4*P-R[I] ; r 3+ r 4*C- r2* P-I [I] 

86: rl-r2*C-r4*P-R[ I+rG] ; r3-r4*C+r 2*P-I [I+r0] 

87: next I ; next J;asp R;next M 

88: for 1 = 0 tc 2*N-1; cll ' 'BI ' ( I , J f N) 

83: if I-U>0;gto "BE" 

SO: it I=J ;g to "INC" 

51: R[I+1J /V(2~N)-P;L[I+1]/y( 2~N)-Z 
92: R[J+1] -R[I + 1] ;I [J+l] +1(1 + 1] 

S3: P-R[J + 1J ; Z-I [J+l] 

94: "INC" :R[I+1] / V( 2"N) -R[I+1] ; 1(1 + 1] /y (2~N) -I ( 1 + 1] 
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95: "EB":next I 

S6: cegjret 

57: "BI" : 0-p2; pl*>p4 

98: for Z=1 tc p3 

59: p4/ 2-*-p4 ; 2*p2-*-p2 

100: if trc ( p4 ) #0 ; p2+l-*-p2 

101: int(p4)-*>p4 

102: next Z;ret 

103: "ABEL " : ina R;lcf 10 ,1$ ,R [ * ] ; R [ 1 ] *N 

104: 1.4*R[1] -1. 8*R[2]-h.4*R[ 3 J -R( 1] 

105: for 1=2 tc 64 
10b: R[I]-M 

107: ,4*N+.2*M-.6*R[I+1] -R[I] ;M-N;next I 

108: for 1=1 tc 64 

109 : R[I ] / (2*V ( (I+i.l) ~ 2-1 * I. ) ) *R [ I ] 

110: for J=I+1 tc 64 

111: R[I]+R[J]/V((J + .l) ■'2-I*L)-R[ I] 

112: next J 

313: if R[L] < .01; 0 — R ( I J 

114: R [ I J / tt -►R [ I J ; asp I;next I 

115: for 1 = 65 tc 512; G-*-R [I ] ; next I; ret 

116: "PLOT" : laf p2,I$,R[*] 

117: 0-*-A;pl-*-B;min(R[*] ) *C ; ma x (R [ * ] ) -»-D 

116: scl A-.l(B-A) ,B-h.G5|E-A) ,C-.1(C-C) ,C+.05(D-C) 

119: B-*-E ; 10-*F ; if flg0;prt A,B,C,D 

120: if pl = 51 2; 6 4*G 

121: if pi = 256; 32 -►G 

122 : if pl= 64; 8-G 

123: ii pl=32;4-*-G 

12 4 : pit B , C , 1 

125: for L=E to 0 ty -G 

126: pit I (E-A) /E+A ,C ,2 

127: pit I (B-A) /E+A ,C+ (D.-C)/ 150 ,2 

128: pit I (B-A) /E+A, C, 2 

129: next I 

130: for 1=0 to F 

131: pit A, I ( B-C) /F+C , 2 

132: pit A+( B-A) /150 , I (D-C)/F+C , 2 

133: pit A, I(D-C)/F+C,2 

134: next I;pen 

135: csiz 1.2,1, .7 

136 : f xq 1 

137: for I=F to 0 by -1 

138: pit A- . 0 7 5 ( E-A ) , I ( D -C ) / F +C , 1 

135: lbl I/F 

140: next I 

141: fxa 0 

142: if p3=2; "LINES/MICRGRADIAN" ♦A$ ; 1/ ( 2*L) ♦L 
143: for 1=0 to E by G 
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144: pit At(I/E-.62cO ( o - A) ,C _ .C25(D — C ) , 1 
145: lbl I* L 

l4o : next I ; it p3 #2; "MICKORADI ANS " -A$ 

147: it p2=2 ; S-*L 

148: pit . 4 ( B-A) +A . C 5 (D -C) +C ; lbl A$ 

149: pit .-.07(B-A)+A,..3(D-C)+C;csiz 1.2,1/. 7, 90 
150: "NORMALIZED IN.TENSI.TY" -*-A$; lbl A$ 

151: csiz . 5 , 1 , 1 . 5 , 0 ; 0*1 

152: pit I , ir.a x ( R ( * ] ) ,1 

153: for 1=1 to pl;plt 1-1 ,R [ I J ; nex t I 

154 : csiz 1 . 2 , 1 * . 7 , 0 

155: pit . 6£ , . SD , 1 ; lol "PLOT OF" 

156: pit .6B,.87D,l;lol 1$ 

157: pit . 6B , . 84D , 1; lbl E$," LASER" 

158: "RET" : pen ;ctg ;ret 0-*Y 
* 175 24 
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